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Abstract 

Generalizing a decision problem for bipartite perfect matching , J. Edmonds introduced 
in |15| the problem (now known as the Edmonds Problem) of deciding if a given linear sub- 
space of M(N) contains a nonsingular matrix, where M(N) stands for the linear space of 
complex N x N matrices. This problem led to many fundamental developments in matroid 
theory etc. 

Classical matching theory can be defined in terms of matrices with nonnegative entries. The 
notion of Positive operator, central in Quantum Theory, is a natural generalization of ma- 
trices with nonnegative entries. (Here operator refers to maps from matrices to matrices.) 
First, we reformulate the Edmonds Problem in terms of of completely positive operators, 
or cquivalently, in terms of bipartite density matrices . It turns out that one of the most 
important cases when Edmonds' problem can be solved in polynomial deterministic time, 
i.e. an intersection of two geometric matroids, corresponds to unentangled (aka separable ) 
bipartite density matrices . We introduce a very general class (or promise ) of linear sub- 
spaces of M(N) on which there exists a polynomial deterministic time algorithm to solve 
Edmonds' problem . 

The algorithm is a thoroughgoing generalization of algorithms in , [SB] , and its analysis 
benefits from an operator analog of permanents, so called Quantum Permanents . Finally, 
we prove that the weak membership problem for the convex set of separable normalized 
bipartite density matrices is NP-HARD. 



1 Introduction and Main Definitions 

Let M(N) be the linear space of N x N complex matrices . The following fundamental problem 
has been posed by J. Edmonds in [T5] : 

Problem 1.1: Given a linear subspace V C M(N) to decide if there exists a nonsingular 
matrix A G V . I 



We will assume throughout the paper that the subspace V is presented as a finite spanning 
A;-tuple of rational matrices S(V) = {A\, A^\{k < N 2 ), i.e . the linear space generated by 
them is equal to V. As usual, the complexity parameter of the input < S(V) > is equal to ( N 
+ "number of bits of entries of matrices Ai, 1 < i < k" ). 

Thus Edmonds' problem is equivalent to checking if the following determinantal polynomial 

P A (xi,...,x k ) = det( ^ XiAi) 

Ki<k 
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is not identically equal to zero. 

This determinantal polynomial can be efficiently evaluated, hence randomized poly-time algo- 
rithms, based on Schwartz's lemma or its recent improvements, are readily available (notice 
that our problem is defined over infinite field with infinite characteristic). 
But for general linear subspaces of M(N), i.e. without an extra assumption (promise), poly-time 
deterministic algorithms are not known and the problem is believed to be "HARD" . 
Like any other homogeneous polynomial, Pa(%1> ■■■■> x k) 1S a weighted sum of monomials of 
degree N, i.e. 

Pa(%1, %k) = a r\,...,r k X\ x 2 ••■ X k i (X) 

(,n,...,r k )eik,N 

where I^jv stands for a set of vectors r = (ri, ■■■,rk) with nonnegative integer components and 
Ei<i<fc1 = N . 

We will make substantial use of the following (Hilbert) norm of determinantal polynomial P(.) 

\\P\\g=- E ki,..,rj 2 r-i!r 2 !...r ifc ! (2) 

(ri,—,r k )eIk,N 

It is easy to show that the determinantal polynomial Pa(xi, %k) = iff -PaC^I) •••> r k) = 
for all (ri, rk) G Ik,N, which amounts to |ifc,jv| = ^N\(k-i)\ com P u t a ti° ns °f determinants . 
We will show that ||-P||q can be evaluated in 0(2^^!) computations of determinants . 
More importantly, ||P||q serves as a natural tool to analyze our main algorithm . 
The algorithm to solve Edmonds' problem, which we introduce and analyze later in the pa- 
per, is a rather thoroughgoing generalization of the recent algorithms |29| . |38| for deciding 
the existence of perfect matchings. They are based on so-called Sinkhorn's iterative scaling. 
The algorithm in [3S] is a greedy version of Sinkhorn's scaling and has been analyzed using 
KLD-divergence; the algorithm in |29| is a standard Sinkhorn's scaling and a "potential" used 
for its analysis is the permanent. Our analysis is a sort of combination of techniques from |29| . 
|38| . Most importantly, ||-P||g can be viewed as a generalization of the permanent. 
The organization of this paper proceeds as follows. In Section 2 we will recall fundamental 
notions from Quantum Information Theory such as bipartite density matrix, positive and com- 
pletely positive operator, separability and entanglement. After that we will rephrase Edmonds' 
problem using those notions and reformulate the famous Edmonds-Rado theorem on the rank 
of intersection of two geometric matroids in terms of the rank non-decreasing property of the 
corresponding (separable) completely positive operator. We will end Section 2 by introducing a 
property, called the Edmonds-Rado property, of linear subspaces of M(N) which allows a poly- 
time deterministic algorithm to solve Edmonds' problem and will explain how is this property 
is related to quantum entanglement. 

In Section 3 we will express G-norm of a determinantal polynomial Pa(%i, i n terms of 

the associated bipartite density matrix, and we will prove various inequalities and properties of 

G-norm which will be needed later on for the analysis of the main algorithm. 

In Section 4 we will introduce and analyze the main algorithm of the paper, Operator Sinkhorn 

Scaling. 

In Section 5 we will apply this algorithm to solve Edmonds' problem for linear subspaces of 
M{N) having the Edmonds-Rado property. In Section 6 we will prove NP-HARDNESS of the 
weak membership problem for the compact convex set of separable normalized density matrices. 
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Finally, in the Conclusion section we will pose several open problems and directions for future 
research. 

We would like to stress that our paper does not contain explicit connections to Quantum Com- 
puting. It rather aims to study quantum entanglement from the point of view of classical 
computational complexity and computational geometry and to use some ideas and structures 
from Quantum Information Theory to construct and analyse classical algorithms. 
The main algorithm of this paper is a third "generation" of scalings applications to computer 
science problems , starting with ( [22] , [HE] ; applied to bipartite perfect matchings and an 
approximation of the permanent ) and f |35| , jSH] ; applied to an approximation of the mixed 
discriminant and mixed volume ) . 

And here it is used to solve very non-trivial , important and seemingly different problem . 



2 Bipartite density matrices, completely positive operators and 
Edmonds Problem 

Definition 2.1: A positive semidefinite matrix pa,b '■ C N <8> C N — ► C N ® C N is called a 
bipartite unnormalized density matrix 

(BUDM ). If tr{pA,B) = 1 then this pa,b is called a bipartite density matrix. 

It is convenient to represent a bipartite pa,b = p(ii,i2,jxi32) as the following block matrix : 



PA,B 



( A x ,x Ai, 2 ... A hN \ 

^2,1 ^.2,2 • • • A.2,N 
\ 4^,1 An,2 ■ ■ ■ A^ y N j 



(3) 



where A hJl =: {p(iuh,h,h) '■ 1 < *2,j2 < N},1 < i\,j x < N. 
A (BUDM ) p is called separable if 

P = P(X,Y) =: x i x i®ViVi> ( 4 ) 

l<i<K 

and entangled otherwise. 

If the vectors Xi, yf, 1 < i < K in (6) are real then p is is called real separable . 
The quantum marginals are defined as pA = J2i<i<N an d 
p B (i,j) =tr(A id );l<i,j <N. 

Next we define the (BUDM ) pa associated with the fc-tuple 

A = (Ax, A k ): 

PA{h,i2,jx, h) =■ A i(h,i2)Ai(ji,j 2 ), (5) 

l<l<k 

where for a complex number z = x + iy its conjugate z = x — iy. 
Rewriting expression (5) in terms of blocks of pa as in (3), we get that 



A i,j = J2 A l^ e ] A l l <^j< N - 



KKk 
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(In quantum physics language, one can view a tuple A = (Ai, ...,Ak) of complex matrices as 
a tuple of unnormalized bipartite "wave functions" ; and (BUDM ) pa as a corresponding 
mixed bipartite state. ) 

We will call (BUDM ) p weakly separable if there exists a separable p', x Y \ with the same 
image as p: Im{p) = Im(p', XY A. 

( Recall that in this finite dimensional case Im(p) is the linear subspace formed by all linear 
combinations of columns of matrix p.) 

A linear operator T : M(N) -> M(N) is called positive if T(X) y for all X y 0, and strictly 
positive if T(X) y atr(X)I for all X y and some a > 0. A positive operator T is called 
completely positive if 

T(X)= Y, AiXA\-Ai,X GM(N) (6) 

l<i<N 2 

Choi's representation of the linear operator T : M(N) — ► M(N) is a block matrix CH(T)ij =: 
T(ejet). The dual to T with respect to the inner product < X, Y >= tr(XY^) is denoted as T*. 
A very useful and easy result of Choi states that T is completely positive iff CH(T) is (BUDM ) . 
Using this natural (linear) correspondence between completely positive operators and (BUDM 
), we will freely "transfer" properties of (BUDM ) to completely positive operators. For 
example, a linear operator T is called separable iff CH(T) is separable, i.e. 

T(Z)=T (XX) (Z)= £ Xl y\Zy lX \ (7) 

l<i<K 

Notice that CH(T (x ,y)) = P(y,x) and T *x,Y) = T ( Y < X )' 

(The components of the vector y are the complex conjugates of corresponding components of 

y)- I 

Remark 2.2: There is a natural (column by column ) correspondence between M(N) and 
C n2 ^C n ®C n . It works as follows 

{A(i,j),l <i,j< N}€M(N)^ 
(A(l, 1), A(l, AO; ; A(l, A), A(N, N)) T e C N * 

In light of definition (2.1), we will represent a linear subspace V C M(N) = C N (g> C N in 
Edmonds Problem as the image of the (BUDM ) p. And as the complexity measure we will 
use the number of bits of (rational) entries of p plus the dimension A. | 

Definition 2.3: A positive linear operator T : M(N) — ► M(N) is called rank non-decreasing 
iff 

Rank{T(X)) > Rank(X) if X y 0; (8) 
and is called indecomposable iff 

Rank(T(X)) > Rank(X) ifX y andl < Rank(X) < N. (9) 

A positive linear operator T : M(N) — ► M{N) is called doubly stochastic iff T(I) = I and 
T*(I) = I ; called e - doubly stochastic iff DS(T) =: tr({T(I) - I) 2 ) + tr((T*(I) - I) 2 ) < e 2 . I 
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The next Proposition(2.4) is a slight generalization of the corresponding result in |29| . 

Proposition 2.4: Doubly stochastic operators are rank non- decreasing. If either T(I) = I or 
T*(I) = I and DS(T) < N' 1 then T is rank non- decreasing. If DS(T) < (2N + l)" 1 then T 
is rank non- decreasing. 

Let us consider a completely positive operator Ta : M{N) — ► M(N), 

T{X) = J2i<i<k AiXA\, and let L(Ai,A 2 , ...,Ak) be a linear subspace of M{N) generated by 
matrices {A^, 1 < i < k}. It is easy to see that if A £ L(Ai, A 2 , A/-) then A(Im(X)) C 
Im(T(X)) for all X y 0. 

Therefore, if L(A%, A2, Ak) contains a nonsingular matrix then the operator T is rank non- 
decreasing. 

This simple observation suggested the following property of linear subspaces of M(N) : 
Edmonds-Rado Property (ERP) 

A linear subspace V = L(A\,A2, A^) has the (ERP) property if the existence of nonsingu- 
lar matrix in V is equivalent to the fact that the associated completely positive operator Ta is 
rank non-decreasing. In other words, a linear subspace V C M(N) has the (ERP) property if 
the fact that all matrices in V are singular is equivalent to the existence of two linear subspaces 
X,Y C C N such dim(Y) < dim(X) and A(X) C Y for all matrices A £ V. 
The main "constructive" result of this paper is that for linear subspaces of M(N) having the 
ERP there is a deterministic poly-time algorithm to solve Edmonds' problem. 
In the rest of this section we will explain why we chose to call this property Edmonds-Rado, 
will describe a rather wide class of linear subspaces with (ERP) property and will give an 
example of a subspace without it. 



2.1 Examples of linear subspaces of M(N) having Edmonds-Rado Property 

Let us first list some obvious but useful facts about the Edmonds-Rado property. 



1. Suppose that V = L{A x ,A 2 ,...,A k ) C M(N) has the (ERP) and C, D £ M(N) are two 
nonsingular matrices. Then linear subspace Vc,d ='■ L(CA\D, CA%D, CA^D) also has 
the (ERP) . 

2. If V = L(A 1 ,A 2 ,...,A k ) C M(N) has the (ERP) then both W =: L(A\, A j 2 , A f k ) and 
V T = L(Af, Al, Al) have the (ERP) . 

3. Any linear subspace V = L(Ax, A 2 , A k ) C M(N) with matrices {^4j,l < i < k} being 
positive semidefinite has the (ERP). 

4. Suppose that linear subspaces V = L(A 1 ,A 2 , ...,A k ) C M(N{) and W = L{B X ,B 2 , ...,B k ) C 
M(N 2 ) both have the (ERP). Define the following matrices Q £ M(Ni+N 2 ), 1 < i < k 
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Then the linear subspace L(Ci, C 2 , C fc ) C M(iVi + N 2 ) also has the (ERP) . 

A particular case of this fact is that any linear subspace of M(N) which has a basis 

consisting of upper diagonal matrices has the (ERP). 

5. Any 1-dimensional subspace of M(N) has the (ERP) property. 

The next theorem gives the most interesting example which motivated the name "Edmonds- 
Rado Property " . Let us first recall one of the most fundamental results in matroids theory, 
i.e. the Edmonds-Rado characterization of the rank of the intersection of two geometric ma- 
troids. 

Definition 2.5: The intersection of two geometric matroids MI(X,Y) = {(xi,yi),l <i< K} 
is a finite family of distinct 2-tuples of non-zero iV-dimensional complex vectors, i.e. Xi, y\ G C N . 
The rank of MI(X, Y), denoted by Rank(MI(X,Y) is the largest integer m such that there 
exist 1 < i\ < ... < i m < K with both sets {x^ , Xi m } and {y^, ■■■,Ui m } being linearly 
independent. I 

The Edmonds-Rado theorem ( |25| ) states (in the much more general situation of the intersec- 
tion of any two matroids with a common ground set) that 



It is easy to see that Rank(MI(X,Y)) is the maximum rank achieved in the linear subspace 
L(xiy\,...,x K y* K ) ; and Rank(MI(X,Y)) = N iff L( Xl y\,... ,xkVk) contains a nonsingular 
matrix. 

Theorem 2.6: Suppose that T : M(N) M(N), T(X) = Ei<j<z AXAj, is a completely pos- 
itive weakly separable operator, i.e. there exists a family of rank one matrices {x\y\, ...,xiyj} C 
M(N) such that L(A 1 , ...,A L ) = L(xiy|, ...,xiy\). 

Then the following conditions are equivalent : 
Fact 1 T is rank non- decreasing. 

Fact 2 The rank of intersection of two geometric matroids MI(X,Y) is equal to N. 

Fact 3 The exists a nonsingular matrix A such that Im{AXA^) C Im(T(X)), X y 0. 

Fact 4 The exists a nonsingular matrix A such that the operator T'(X) = T(X) — AXA^ is 
completely positive. 



Rank{MI(X,Y)) = 
mm Sc {i 2,.. .,K}dimL(xi;i G S) + dimL(yj; j G S) 



(10) 
(11) 



Proof: [2 =^ 1] 

Suppose that the rank of MI(X, Y) is equal to N . Then 



RankT(X) = dim(L(xf,i 



G S)) where S =: {i : y]X yi / 0} 
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As dim(L{yj\ j G S) < dim(Ker(X)) = N — Rank(X) hence, from the Edmonds-Rado Theo- 
rem we get that 

Rank{T(X)) > N - (N - Rank(X)) = Rank(X). 



It follows from the Edmonds-Rado Theorem that the rank of MI(X, Y) is equal to N. All 

other "equivalences" follow now directly. 

I 

Remark 2.7: Theorem 2.6 makes the Edmonds-Rado theorem sound like Hall's theorem on 
bipartite perfect matchings. 

Indeed, consider a weighted incidence matrix Ay of a bipartite graph T, i.e.^4r(i, j) > if i 

from the first part is adjacent to j from the second part and equal to zero otherwise. Then 

Hall's theorem can be immediately reformulated as follows : 

A perfect matching, which is just a permutation in this bipartite case, exists iff 

|^4r^|+ > \ x \+ for any vector x with nonnegative entries, where stands for the number of 

positive entries of a vector x. 

All known algorithms (for instance, linear programming based on [25]) to compute the rank 
of the intersection of two geometric matroids require an explicit knowledge of pairs of vectors 
(xi,yi), or, in other words, an explicit representation of the rank one basis {xiyj, 1 < i < I}. 
The algorithm in this paper requires only a promise that such a rank one basis (not necessarily 
rational!) does exist. | 

Another example comes from Consider pairs of matrices (Ai,Bi G M(N); 1 < i < K. 
Let Vi C M{N) be the linear subspace of all matrix solutions of the equation XAi = BiX. 
One of the problems solved in JS] is to decide if W = V\ fl ... D Vk contains a nonsingular 
matrix. 

It is not clear to the author whether the class of such linear subspaces W satisfies the (ERP) 
property. 

But suppose that A\ is similar to B\ (V\ contains a nonsingular matrix ) and, additionally, 
assume that dim(Ker(Ai — XI) = dim(Ker(B\ — XI) < 1 for all complex A G C. 
(I.e. just one Jordan block for each eigenvalue.) 

It is not difficult to show that in this case there exist two nonsingular matrices D, Q and upper 
diagonal matrices (Ui, U r ) such that V\ = L(DUiQ, DU r Q). It follows, using 
Facts (1, 4 ) above, that V\ as well as any of its linear subspaces satisfy (ERP). 



Example 2.8: Consider the following completely positive doubly stochastic operator Sk^ : 
M(3) -> M(3) : 




L{yj\ j G S) L . Then 

dim(L(xi : i G S)) > 
Rank(T(P)) > Rank(P) = N - dim(L( yj ;j G S)). 



Sk 3 (X) = -{A {1>2) XA\ 12) +A {1>3) XA\ li3) +A {2>3) XA\ 23) ) 



(12) 
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Here {Auj\, 1 < i < j < 3} is a standard basis in the linear subspace K(3) C M(3) consisting 
of all skew-symmetric matrices, i.e. A(ij) ='■ ejej — eje\ and {ej,l < i < 3} is a standard 
orthonormal basis in C 3 . 

It is clear that all 3 x 3 skew-symmetric matrices are singular. As Sk^ is a completely positive 
doubly stochastic operator, and, thus, is rank non-decreasing, therefore K(3) C M(3) is an 
example of a linear subspace not having (ERP) property. 
More "exotic" properties of this operator can be found in [Jj. | 

3 Quantum permanents and G-norms of determinantal polyno- 
mials 

Consider a fc-tuple of N x N complex matrices A = (A\, A^). Our first goal here is to 
express the square of the G-norm of a determinantal polynomial Pa( x i, ■■■,xi c ) in terms of the 
associated bipartite density matrix (BUDM ) p&_, which is defined as in (5). 
Consider an A-tuple of complex N x N matrices, B = (Bi, -Bjy). Recall that the mixed 
discriminant M(B) = M(B\, ...,Bn) is defined as follows : 

M(B U ...B N ) = dx ® n QxN det(xi£i + .... + x N B N ). (13) 

Or equivalently : 

N 

M(B 1 ,...B N )= ]T (-lr^n^KCUW), (14) 

where S n is the symmetric group, i.e. the group of all permutations of the set {1, 2, • • • , N}. If 
matrices Bi,l < i < N are diagonal then their mixed discriminant is equal to the corresponding 
permanent (JSHl)- 

Definition 3.1: Let us consider a block matrix p as in (3) (not necessarily positive semidefi- 
nite). We define the quantum permanent, QP(p), by the following equivalent formulas : 

QP(p) =: £ (-ir^M(A lMl) ,...,A NMN) ); (15) 



t1,T2,T3,T4£Sn 

N 

Y[rho(Ti(i),T 2 (i),T 3 (i),T 4 (i)). (16) 
i=i 

I 
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Straight from this definition, we get the following inner product formula for quantum perma- 
nents : 

QP(p)=<p® N Z,Z>, (17) 

where p® N stands for a tensor product of N copies of p, < ., . > is a standard inner product 
and 

if jW = T k (i)(l <i<N)-T k £ S N (k = 1,2) and zero otherwise. 



Remark 3.2: Notice that the equality (17) implies that if 
Pi h P2 h then QP(pi) > QP(p 2 ) > . 

The standard norm of ./V 2Ar -dimensional vector Z defined above is equal to 1. Thus, if p is a 
normalized bipartite density matrix then QP{p) can be viewed as the probability of a particular 
outcome of some (von Neumann) measurement. Unfortunately, in this case QP{p) < jpj-. 
Consider an arbitrary permutation a £ S4 and for a block matrix (or tensor ) p = {p(h ,i2,hiU)'A < 
h,i2,i3,U < N} define p a = {p(i a (i), V( 2 ), V(3)> V(4)}- It is easy to see that QP(p) = QP{p a ). 
Another simple but important fact about quantum permanents is the following identity : 

QPdAi A 2 )p(A 3 A A )) = detiA^AsA^QPip) (18) 

The author clearly (and sympathetically ) realizes that some readers might object to (or ridicule) 
the name "quantum permanent". The next example, hopefully, will explain possible motiva- 
tions. I 

Example 3.3: Let us present a few cases when Quantum Permanents can be computed "exactly 
" . They will also illustrate how universal this new notion is. 

1. Let p A ,B be a product state, i.e. p A ,B = C®D. Then QP(C D) = N\Det(C)Det(D). 

2. Let pa,b be a pure state, i.e. there exists a matrix (R = R(i,j) : 1 < i,j < N) 
such that PA,B(k,i2,ji,j2) = R(ii,i2)R(ji,j2)- 

In this case QP{pa,b) = N\\Det{R)\ 2 . 

3. Define blocks of pa,b as Aij = R(i,j)eie\. 
Then QP( Pa ,b) = Per{R). 

I 

The following propositions provide important upper bounds for quantum permanents of positive 
semidefinite matrices. 

Proposition 3.4: Suppose that pa,b is a (T3UDM ). Then 

max \M(A 1M1) ,...,A NMN) )\ = 

M{A 1A ,...,A NjN ) (19) 
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Proof: For r, a G Sn define a matrix 

B T . a =: -A t (1),<t(1) ® ^t(2),<t(2) ® ■•• ® ^r(7V),<r(Ar) 

Since /oa.b is positive semidefinite hence the block matrix {-B t ,<t : t, cr G S^r} is also positive 
semidefinite. It is well known ( and easy to prove that 

M(A U .., A N ) = tr{{Ax ® ... ® ^yW) 

for some universal TV^-dimensional vector V . 
It follows that the following Nl x Nl matrix C 

C r ,a = tr{B T!a VV<) = M(A T ( 1 ) j(T ( 1) , ^4 r (2), (7 (2) 1 •••) At(N),<t(N)) 

is also positive semidefinite. Thus 

\C Tt a\ < (C T>T C a , a )3 =M(A ltl ,...,A N>N ) 

I 

Corollary 3.5: If p A ,B is (BUDM ) then 

QP(pa,b) < N\M(A 1>U A NtN ) < N\Det( PA ). (20) 

The permanental part of Example (3. 3) shows that Nl is the exact constant in both parts of (20), 
i.e. if blocks Aij = ejej, 1 < i,j < N 

then QP(pa,b) = Nl and M(A hl , A N}N ) = Det(p A ) = 1. 

The next proposition follows from Hadamard's inequality : 
if X y is N x N matrix then Det(X) < l[f =1 X(i, i). 

Proposition 3.6: If X >- then the following inequality holds : 

K 

DetC^XiyjXyixl) > 

i=l 

Det{X)MP {xy) . (21) 

Corollary 3.7: Suppose that a separable (T3UDM ) p A ,B is Choi's representation of the 

completely positive operator T. 

Then for all X y the following inequality holds : 

Det(T(X)) > QP(p AB )Det(X) (22) 

Since p A = T(I), hence QP(p a ^b) < Det(p A ) in the separable case. 
(Notice that Corollary 3.5 provides an example of an entangled 
fBUDM ) which does not satisfy (22) .) 
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Finally, in the next theorem we connect quantum permanents with G-norms of determinan- 
tal polynomials. 



Theorem 3.8: 



1. Consider an arbitrary polynomial 

P(x 1 ,x 2 ,...,x k ) = a n,...,r k x r 1 1 x r 2 2 ...x r k k ,\F\ < oo 

(ri,...,r fe )eF 

where F is some finite set of vectors with nonnegative integer components and define its 
G -norm as follows 

\\P\\g= : k 1 ,...,rj 2 ri!r 2 !...r fc ! 

{n,...,r k )eF 

Then the following identity holds : 

\\P\\G = E il> ...^(\P^i,...^ k )\ 2 ), (23) 

where (£i, ...,^k) are independent identically distributed zero mean gaussian complex ran- 
dom variables and the covariance matrix o/£i, viewed as a 2-dimensional real vector, is 
equal to \l. 

2. Consider a k-tuple of N x N complex matrices A = (A±, A^) and the corresponding 
determinantal polynomial Pa(xi, x^) =: det(^ 1<i<fc XiAi). Then the following identity 
holds 

\\Pa\\1 = QP(pa) (24) 



Proof: The proof is in Appendix 1. | 



Remark 3.9: Theorem 3.8, more precisely the combinations of its two parts, can be viewed as 
a generalization of the famous Wick formula |28| 

It seems reasonable to predict that formula(24) might be of use in the combinatorics described 
in HE!- 

It is well known (see, for instance, JJj ) that the mixed discriminant M(A\, An) can be 
evaluated by computing 2 N determinants. Therefore there the quantum permanent QP{p) 
can be evaluated by computing N\2 N determinants. Now, formula (24) suggests the following 
algorithm to compute 

l|det(Ei<i< fc ^)llG : 

first, construct the associated bipartite density matrix p&_, which will take 0{N 4L k) additions 
and multiplications ; secondly, compute QP(pa)- 

Total cost is Cost(N) = 0(N12 N iV 3 ). On the other hand, just the number of monomials in 
detQCi<i<fc ^Ai) is equal to \h,N\ = ^N\(k-i)\ • If k ~ 1 = aN2 then I-^mH - • Thus, 

\h n\ a N N 2N (ae 2 ) N 



^ > ; „ „ >f 



Cost(N) ~ 0(Nl 2 2 N N 3 ) ~ 2 N N* 
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We conclude that if a > \ our approach is exponentially faster than the "naive" one, i.e. than 
evaluating det(^i<j<fe x i^i) a t ah vectors (xi, ...,Xk) £ Ik,N- 

Our approach provides an 0(N12 N N 3 ) - step deterministic algorithm to solve a general case of 

Edmonds' Problem . 

I 



4 Operator Sinkhorn's iterative scaling 



Recall that for a square matrix A = {a^ : 1 < i, j < N} row scaling is defined as 

column scaling as C(A) = — } assuming that all denominators are nonzero. 

The iterative process ...CRCR(A) is called Sinkhorn's iterative scaling (SI). There are two 
main, well known, properties of this iterative process, which we will generalize to positive 
operators. 



Proposition 4.1: 



1. Suppose that A = {aij > : 1 < i, j < N}. Then (SI) converges iff A is matching, i.e., 
there exists a permutation 7r such that a ij7r (j) > (1 < i < N). 

2. If A is indecomposable, i.e., A has a doubly-stochastic pattern and is fully indecomposable 
in the usual sense, then (SI) converges exponentially fast. Also in this case there exist 
unique positive diagonal matrices Di,D2,det(D2) = 1 such that the matrix D^^^AD^ 1 is 
doubly stochastic. 

Definition 4.2: [Operator scaling ] Consider a positive linear operator T : M(N) — > M(N). 
Define a new positive operator, Operator scaling, S , c 1 ,c , 2 (^) as : 

S Cl ,c 2 (T)(X) =: dT(ClxC 2 )Cl (25) 

Assuming that both T(I) and T*(I) are nonsingular we define analogs of row and column 
scalings : 

^) = V)-i/ (T) ' c(T) = V(/)-i (T) (26) 

Operator Sinkhorn's iterative scaling (OSI) is the iterative process ...CRCR(T) | 



Remark 4.3: Using Choi's representation of the operator T as in Definitional), we can de- 
fine analogs of operator scaling (which are exactly so called local transformations in Quantum 
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Information Theory) and (OSI) in terms of (BUDM ) : 

Sc\,c 2 (pa,b) = Ci ® C 2 (pa,b)C\ ® C 2 f ; 

_ i _ i 

R(pa,b) = Pa 2 ® I{pa,b)p a 2 ® ^ 

C(pa,b) = I ® PbHpa,b)I ® Pb* ■ (27) 

The standard ("classical") Sinkhorn's iterative scaling is a particular case of Operator Sinkhorn's 
iterative scaling (OSI) when the initial Choi's representation of the operator T is a diagonal 
(BUDM ) . | 

Let us introduce a class of locally scalable functionals (LSF ) defined on a set of positive linear 
operators, i.e. functionals satisfying the following identity : 

<p(S Cl MT)) = DetidC^DetiCiCiMT) (28) 

We will call an (LSF ) bounded if there exists a function / such that \<p(T)\ < f(tr(T(I)). 
It is clear that bounded (LSF ) are natural "potentials" for analyzing (OSI). Indeed, Let 
T n ,To = T be a trajectory of (OSI). T is a positive linear operator. Then T,i{I) = I for odd i 
and T 2i (I)* = I,i>l. Thus if tp(.) is (LSF ) then 

<p(T i+1 ) = a(i)<p{Tj,a(i) = Detail))- 1 if i is odd , 

a(i) = Detail))' 1 if % > is even. (29) 

As tr(Ti(I)) = tr(T*(I)) = N,i > 0, thus by the arithmetic/geometric means inequality we 
have that |<p(Tj + i)| > |<£>(Tj)| and if ip(.) is bounded and |y(T)| / then DS(T n ) converges to 
zero. 

To prove a generalization of Statement 1 in Prop. (4.1) we need to "invent" a bounded (LSF 
) (p{.) such that <p(T) ^ iff the operator T is rank non-decreasing. We call such functionals 
"responsible for matching". It follows from (10) and (20) that QP(CH(T)) is a bounded (LSF 
). Thus if QP(CH(T)) / then DS(T n ) converges to zero and, by Prop. (2.4), T is rank 
non-decreasing. On the other hand, QP(CH(Skz)) = and Sks is rank non-decreasing (even 
indecomposable ). This is another "strangeness" of entangled operators. We wonder if it is 
possible to have a "nice", say polynomial with integer coefficients, responsible for matching 
(LSF ) ? We introduce below a responsible for matching bounded (LSF ) which is continuous 
but non-differentiable. 

Definition 4.4: For a positive operator T : M(N) — > M(N), we define its capacity as 
Cap{T) = mf{Det(T(X)) : X y 0,Det(X) = 1}. I 

It is easy to see that Cap(T) is (LSF ). 
Since Cap(T) < Det(T{I)) < (^ r ™ 1 ) JV , 
hence Cap(T) is a bounded (LSF ). 
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Lemma 4.5: A positive operator T : M(N) — ► M(N) is rank non- decreasing iff Cap(T) > 0. 

Proof: Let us fix an orthonormal basis (unitary matrix) U = {u\, m C N and associate 

with a positive operator T the following positive operator : 

T V (X)=: £ T(u iU l)tr(X Ui ul). (30) 

l<i<N 

(In physics terms, Tjj represents decoherence with respect to the basis U, i.e. in this basis 
applying Tjj to matrix X is the same as applying T to the diagonal restriction of X. ) 
It is easy to see that a positive operator T is rank non-decreasing iff the operators Tjj are rank 
non-decreasing for all unitary U. 

And for fixed U all properties of Tjj are defined by the following Af-tuple o£ N x N positive 
semidefinite matrices : 

A T ,u =■ (T(uiu\), ...,T(u N vJ N ). (31) 

Importantly for us, Tjj is rank non-decreasing iff the mixed discriminant M(T(uiu\), T(ujv^jv)) ^ 
0. 

Define the capacity of At,u, 

Cap(A TtU ) =: 
\ni{Det(Y.i<i<N T ( u i u] i)li) : H > )Ili<i<7v7i = li- 
lt is clear from the definitions that Cap{T) is equal to infimum of Cap(AT,u) ° ver an unitary 
U. 

One of the main results of |35| states that 

M(At.u) =■ M(T(u lU \),...,T{uW N ))<Cap(A TiU )< 
N N 

<j^M(T(u lU \),...,T(uW N ))- (32) 

As the mixed discriminant is a continuous (analytic) functional and the group SU(N) of unitary 
matrices is compact, we get the next inequality: 

min M(Atu) < Cap(T) < — - min M(A TU ) (33) 

The last inequality proves that Cap(T) > iff positive operatorT is rank non-decreasing. | 

So, the capacity is a bounded (LSF ) responsible for matching, which proves the next 
theorem : 

Theorem 4.6: 

1. Let T h ,Tq = T be a trajectory of (OSI), where T is a positive linear operator. Then 
DS(T n ) converges to zero iffT is rank non- decreasing. 

2. A positive linear operator T : M(N) — > M(N) is rank non- decreasing iff for all e > 
there exists an e-doubly stochastic operator scaling of T . 
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3. A positive linear operator T is rank non- decreasing iff there exists i -doubly stochastic 
operator scaling of T. 

The next theorem generalizes second part of Prop. (4.1) and is proved on almost the same lines 
as Lemmas 24,25,26,27 in 

Theorem 4.7: 

1. There exist nonsingular matrices C\,C% such that S^c^T) is doubly stochastic iff the 
infimum in Definition 4-4 * s attained. 

Moreover, if Cap(T) = Det{T(C)) where C y 0,Det(C) = 1 
then S -i i (T) is doubly stochastic. 

Positive operator T is indecomposable iff the infimum in Definition 4-4 ^ s attained and 
unique. 

2. A doubly stochastic operator T is indecomposable iff 

tr(T(X)) 2 < a tr(X) 2 for some < a < 1 and all traceless hermitian matrices X. 

3. If a positive operator T is indecomposable then DS(T n ) converges to zero with the expo- 
nential rate, i.e. DS(T n ) < Ka n for some K and < a < 1. 

Remark 4.8: Consider an N x TV matrix A with nonnegative entries. Similarly to (30), define 
its capacity as follows : 

Cap(A) = inf{ JJ (Ax) l : x { > 0, 1 < i < N; J[ x l = 1} 

l<i<N l<i<N 

Recall that the KLD-divergence between two matrices is defined as 

KLD(A\\B)= J2 S(^j)log(^M) 
l<i,j<N A ^3> 

It is easy to prove (see, for instance, [SH] ) that 

-\og{Cap{A)) = mi{KLD(A\\B) : B e D N ], 

where D^ is the convex compact set of N x N doubly stochastic matrices. 
Of course, there is a quantum analog of KLD-divergence, the so called von Neumann divergence. 
It is not clear whether there exists a similar "quantum" characterization of the capacity of 
completely positive operators. 

The inequality (20) can be strengthen to the following one : 

QP{CH{T)) < N\Cap{T) 
And AM is also an exact constant in this inequality above. I 
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5 Polynomial time deterministic algorithm for the Edmonds 
Problem 



Let us consider the following three properties of (BUDM ) pab- ( We will view this pab a s 
Choi's representation of a completely positive operator T, i.e. pa,b = CH(T). ) 

PI Im{pA,B) contains a nonsingular matrix. 
P2 The Quantum permanent QP(pa,b) > 0. 
P3 Operator T is rank non-decreasing. 



Part 2 of theorem (3.8) proves that PI 44> P2 and Example(2.8) illustrated that that the impli- 
cation P2 => P3 is strict. It is not clear whether either PI or P3 can be checked in deterministic 
polynomial time. 

Next, we will describe and analyze Polynomial time deterministic algorithm to check whether 
P3 holds provided that it is promised that Im{pA,B) , viewed as a linear subspace of M(N), 
has the Edmonds-Rado Property . Or, in other words, that it is promised that PI P3. 
In terms of Operator Sinkhorn's iterative scaling (OSI) we need to check if there exists n 
such that DS(T n ) < If L =: min{n : DS(T n ) < -^} is bounded by a polynomial in N 
and number of bits of pa,b then we have a Polynomial time Deterministic algorithm to solve 
Edmonds' problem provided that it is promised that PI 4=> P3. Algorithms of this kind for 
"classical" matching problem appeared independently in and jSH]- Iri the "classical" case 
they are just another, conceptually simple, but far from optimal, poly-time algorithms to check 
whether a perfect matching exists. But in this general Edmonds Problem setting, our, Operator 
Sinkhorn's iterative scaling based approach seems perhaps to be the only possibility. 
Assume, without loss of generality, that all entries of pa,b are integer numbers and their max- 
imum magnitude is M. Then Det(pA) < (MN) N by Hadamard's inequality. If QP(pa,b) > 
then necessary QP(pa,b) > 1 for it is an integer number. Thus 

QP(CH( Tl )) = QPVfF)) > N 
Det(pA) 



Each nth iteration (n < L ) after the first one will multiply the Quantum permanent 
by Det(Xy l , where X >- 0,tr(X) = N and tr((X - I) 2 ) > ^. Using results from [15], 
Det(X)- 1 > (1 - 3^) _1 =: 6. Putting all this together, we get the following upper bound on 
L, the number of steps in (OSI) to reach the "boundary" DS(T n ) < jt : 



sL ^ QP(CH( T L )) 
(MN)" 



5 l < (34 ) 



It follows from (20) that QP(CH(T L )) < N\ 



Taking logarithms we get that 

L 3N{Nln(N) + N(ln(N) + ln(M)); (35) 
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Thus L is polynomial in the dimension N and the number of bits log(M). 
To finish our analysis, we need to evaluate the complexity of each step of (OSI). 
Recall that T n (X) = L n (T(RlXR n ))Ll 
for some nonsingular matrices L n and R n , 
T n {I) = L n (T(RiR n ))Ll and T*{I) = R n (T*(LtL n ))Rl 
To evaluate DS(T n ) we need to compute tr((T*(I) — I) 2 ) for odd 
n and tr((T n (I) — I) 2 ) for even n. 

Define P n = L^Ln, Q n = R^Rn- It is easy to see that the matrix T n (I) is similar to P n T(Q n ), 
and T*(I) is similar to Q n T*(P n ). 

As traces of similar matrices are equal, to evaluate DS(T n ) it is sufficient to compute matrices 

Pn ? Qn ■ 

But, P n+l = (T{Q n ))~ l and Q n+1 = (T*{P n ))- x . 

And this leads to standard rational matrix operations with 0(N 3 ) per iteration in (OSI). 
Notice that our original definition of (OSI) requires computation of an operator square root. 
It can be replaced by the Cholesky factorization, which still requires computing scalar square 
roots. But our final algorithm is rational! 

Remark 5.1: To ensure that all the matrices we need to invert along the algorithm (OSI) are 
nonsingular indeed it is sufficient that both T(I) y and T*(I) y (strictly positive definite) . 
It is easy to see that if positive operator T is rank non-decreasing then its dual T* is also rank 
non-decreasing . Thus if positive operator T is rank non-decreasing then necessarily T(X) y 
and T*{X) y for all X y . I 

6 Weak Membership Problem for the convex compact set of 
normalized bipartite separable density matrices is NP-HARD 

One of the main research activities in Quantum Information Theory is a search for " operational" 
criterion for the separability. We will show in this section that, in a sense defined below, the 
problem is NP-HARD even for bipartite normalized density matrices provided that each part 
is large (each "particle" has large number of levels). First, we need to recall some basic notions 
from computational convex geometry. 

6.1 Algorithmic aspects of convex sets 

We will follow 

Definition 6.1: A proper (i.e. with nonempty interior) convex set K C R n is called well- 
bounded a-centered if there exist a rational vector a £ K and positive (rational) numbers r, R 
such that B(a,r) C K and K C B(a,R) (here B(a,r) = {x : \\x — a\\ < r} and |.| is a standard 
euclidean norm in R n ) . The encoding length of such a convex set K is 

< K >= n+ < r > + < R > + < a > , 
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where <r>,<R>,<a> are the number of bits of corresponding rational numbers and 
rational vector. 

Following [25] we define S(K, 5) as a union of all <5-balls with centers belonging to K ; and 

S(K, -5) = {x G K : B(x, 5) C K}. 

I 

Definition 6.2: The Weak Membership Problem (WMEM(K,y,5)) is defined as follows : 
Given a rational vector y G i? n and a rational number 5 > either 

(i) assert that y G S (K, 5) , or 

(ii) assert that y S(K, —5). 

The Weak Validity Problem (WVAL(K, c, 7, 5)) is defined as follows : 
Given a rational vector c G rational number 7 and a rational number 5 > either 

(i) assert that < c, x >=: c T x < 7 + <5 for all x G ^(i^, —5), or 

(ii) assert that c T x > 7 — (5 for some x G 5). 
I 

Remark 6.3: Define M(K,c) =: max xG ^ < c, x > . It is easy to see that 

M(K,c)>M(S(K,-S),c) >M(K,c)-\\c\\Sf, 
M(K,c) < M(S(K,5),c) < M(K,c) + ||c||<5 

I 

Recall that the seminal Yudin-Nemirovski theorem ([Hj, implies that if there exists a 

deterministic algorithm solving WMEM(K, y, 5) in Poly{< K> + <y> + <5>) steps then 
there exists a deterministic algorithm solving WVAL(K, c, 7, 5) in Poly{< K> + <c> + < 
5 > + < 7 >) steps. 

Let us denote as SEP(M, N) a compact convex set of separable density matrices pa,b '■ C M ® 
C iv ^ C Af c ,jv ) tr( P(A)B ) = 1, M > N. Recall that 

SEP(M,N) = 
CO({xxt ®yyt . 3. G c M ,i/ G C^; ||x|| = ||y|| = 1}), 

where CO(X) stands for the convex hull generated by a set X. 

Our goal is to prove that the Weak Membership Problem for SEP(M, N) is NP-HARD. As we 
are going to use the Yudin-Nemirovski theorem, it is sufficient to prove that W VAL(SEP(M, N),c, 7, 5) 
is NP-HARD with respect to the complexity measure (M+ <c> + <<5> + <7>) and to 
show that < SEP(M, N) > is polynomial in M. 

6.2 Geometry of SEP(M, N) 

First, SEP(M, N) can be viewed as a compact convex subset of the hyperplane in R D ,D =: 
N 2 M 2 . The standard euclidean norm in R N M corresponds to the Frobenius norm for den- 
sity matrices, i.e. ||p||f = t r (pp^)- The matrix j^jl G SEP(M, N) and Htv^j^ — xx^ (g> 
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VU^Wf = Y "TT" ^ ^ f° r an norm one vectors x,y. Thus SEP(M, N) is covered by the ball 
The following result was recently proved in |4U| . 

Theorem 6.4: Let A be a block hermitian matrix as in (5). Iftr(A) = and [| A||i? < rj(D-i ) 
then the the block matrix j^I + A is separable. 

Summarizing, we get that for D = MN 

(balls are restricted to the corresponding hyperplane ) and conclude that < SEP(M,N) >< 
Poly(MN). It is left to prove that WVAL(SEP(M, N), c, 7, 5) is NP-HARD with respect to 
the complexity measure (MN+ <c> + <5> + <7>). 



6.3 Proof of Hardness 



Let us consider the following hermitian block matrix : 



/ 



C 





Ax 








A M -i \ 


J 



(36) 



i.e. its (i,j) blocks are zero if either i 7^ 1 or j 7^ 1 and (1, 1) block is also zero ; Ai, ...,Am-i 
are real symmetric N x N matrices. 



Proposition 6.5: 

max pe5jB p (Mi7V) (tr(Cp)) 2 = 
max J/gjR iv i ||j / || =1 Ei<i<A'/-i(y T - 4 iy) 2 - 



Proof: First, by linearity and the fact that the set of extreme points 

Ext(SEP(M,N)) = 
{xx^ (g> yy^ : x £ C M , y G C N ; \\x\\ = \\y\\ = 1} 

we get that 

max p6SE? ( M , A f) tr(Cp) = 

™® x -xxt®yyi-.x€C M ,y€C N -,\\x\\ = \\y\\=l tr{C(xX^ (g> yy*)). 
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But tr(C(yy^ (8> xx')) = tr(A(y)xx^), where real symmetric M x M matrix A(y) is defined as 
follows : 

/ ax ... OAf-i \ 
ai ... 



Mv) 







(H = tr(Aiyyt), 1 < i < M - 1. 
Thus 



max peSi? p(M,jv) ir(Cp) = 
max w t® :Ea; t :a; eG^,2/eC^;[| a; ||=|| ? /||=i ^(CCsxt <g) yyt)) = 
maxi^i^! A ma:c A(y). 

(Above A maxA{y) is the maximum eigenvalue of A(y)) 

It is easy to see A(y) has only two real non-zero eigenvalues 

(d, -d), where d = T,i<i<M-i( tr ( A iyy' ( )) 2 ■ 

As Ai,l < i < N — 1 are real symmetric matrices we finally get that 

max (tr(Cp)) 2 = max (y T A;y) 2 . 



Proposition(6.5) and Remark(6.3) suggest that in order to prove NP-HARDness of 
WVAL(SEP(M, N),c, 7, S) with respect to the complexity measure M+ <c> + <5> + < 
7 > it is sufficient to prove that the following problem is NP-HARD : 



Definition 6.6: (RSDF problem) Given k I xl real rational symmetric matrices (Ai, 1 < i < I) 
and rational numbers (7, 5) to check whether 

7 + 5 > max xeR i M=1 f(x) > 7 - 6, f{x) =: ^ (x T Aixf. 

l<i<l 

respect to the complexity measure 

(Ik + Ei<i<i < Ai > + < 8 > + < 7 >). I 

It was shown in [T3j. bv a reduction from KNAPSACK, that the RSDF problem is NP-HARD 
provided 

We summarize all this in the following theorem 

Theorem 6.7: The Weak Membership Problem for SEP(M, N) is NP-HARD if N < M < 

m^n + 2 . 

Remark 6.8: It is easy exercise to prove that (BUDM ) pa,b written in block form (3) is 
real separable iff it is separable and all the blocks in (3) are real symmetric matrices. It follows 
that, with obvious modifications, Theorem 6.7 is valid for real separability too. 
The construction (37) was inspired by Arkadi Nemirovski's proof of the NP-HARDness of 
checking the positivity of a given operator I 
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7 Concluding Remarks 



Many ideas of this paper were suggested by [33] . The world of mathematical interconnections 
is very unpredictable (and thus is so exciting). The main technical result in a very recent 
breakthrough in Communication Complexity 39 is a rediscovery of particular, rank one, case 
of a general, matrix tuples scaling, result proved in |35| with much simpler proof than in |39j . 
Perhaps this paper will produce something new in Quantum Communication Complexity. 
We still don't know whether there is a deterministic poly-time algorithm to check whether a 
given completely positive operator is rank non-decreasing. And this question is related to lower 
bounds on Cap(T) provided that Choi's representation CH(T) is an integer semidefinite ma- 
trix. 

Theorem(6.7) together with other results from our paper gives a new, classical complexity 
based, insight on the nature of quantum entanglement and, in a sense, closes a long line of 
research in Quantum Information Theory. 

Also, this paper suggests a new way to look at "the worst entangled" bipartite density matrices 
(or completely positive operators). For instance, the operator Sk^ from Example (2.8) seems 
to be "the worst entangled" and it is not surprising that it appears in many counterexamples. 
We hope that the constructions introduced in this paper, especially the Quantum Permanent, 
will have a promising future. 

We think, that in general, mixed discriminants and mixed volumes should be studied and 
used more enthusiastically in the Quantum context. After all, they are noncommutative gen- 
eralizations of the permanent.... 

The G-norm defined in (2) appears in this paper mainly because of formula (24). It is called 
by some authors ( |2j ) Bombieri's norm (see also, [2], [I], |SJ )• 

Also, the G-norm arises naturally in quantum optics and the study of quantum harmonic os- 
cillators. 

This norm satisfies some remarkable properties ( [3J, [1]) which, we think, can be used in quan- 
tum/linear optics computing research. 

Combining formulas (23) and (24) , one gets an unbiased nonnegative valued random estimator 

for quantum permanents of bipartite unnormalized density matrices. But, as indicated in 

it behaves rather badly for the entangled bipartite unnormalized density matrices. 

From the other hand, there is a hope, pending on a proof of a "third" generation of van der 

Waerden 

conjecture (([23]) [22], I2B )' (I3S], [2H])) > to have even a deterministic polynomial time al- 
gorithm to approximate within a simply exponential factor quantum permanents of separable 
unnormalized bipartite density matrices (more details on this matter can be found in 7 ). It is 
my great pleasure to thank my LANL colleagues Manny Knill and Howard Barnum . 
Many thanks to Marek Karpinski and Alex Samorodnitsky for their comments on this paper . 
Finally, I would like to thank Arkadi Nemirovski for many enlightening discussions. 

A Proof of Theorem (3.8) and a permanental corollary 

The main goal of this Appendix is a "direct proof of formula (24) . A much shorter probabilistic 
proof is presented in Appendix C . 
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Proof: [ Proof of formula (23) ] 

It is sufficent to prove that for any monomial 

J ... J \z[\..zl k \ 2 e-( x i +y i\..e-( x l +y lUx 1 dy 1 ...dx k dy k = r 1 !r 2 !...r fe !(z i = xi + iy h l <l<k). 

(37) 

And that distinct monomials are orthogonal , i.e. 

j ... j (^^..^ fc ^...zJ fc )e-^ +w ?)...e-^ + ^)cteidyi...tte fc dy fc = 0(r + h) (38) 

Notice that both 2/c-dimensional integrals (37) and (38) are products of corresponding 2- 
dimensional integrals . Thus (37) is reduced to the fact that 

i j J(xj + ylf^e'^^dx^ = n!. 

Using polar coordinates in a standard way , we get that 

Similarly (38) is reduced to 

/ J(xi + iyi) m {x\ + yl) k e-^ + y^d Xl d yi = 0, 

where m is positive integer and k is nonnegative integer . 
But 



J J(x! + iyi) m {x\ + yl) k e- x ' +y U Xl d yi = R 2k e~ R2 (J** e~ im ^d(j))dR = 0. 



Proof: [Proof of formula (24) ] 

First , let us recall how coefficients of det(J2i<i< k XiAi) can be expressed in terms of the 
corresponding mixed discriminants . Let us associate a vector r G I k ,N an -/V-tuple of N x N 
complex matrices B r consisting of copies of Ai(l < i < k) . 
Notice that 

B r = (B u B N ); B t e {A lt A k }, 1 < i < k. 

It is well known and easy to check that for this particular determinantal polynomial its coeffi- 
cients satisfy the following identities : 

; (n,...,r k ) e I k ,N (39) 



r\\r2\-..r k \ 



We already defined mixed discriminants by two equivalent formulas (13) , (14) . The next 
equivalent definition is handy for our proof : 

M(B 1 ,...B N ) = det([B 1 (e a{1) \B 2 (e a{2) \...\B N (e a(N) ]). (40) 

creS N 
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In the formula( 40) above , (e±, ejv) is a canonical basis in C N , and for a N x TV complex 
matrix 5 a column vector B{ei) is an i-th column of i? . 

We will use in this proof three basic elementary facts about mixed discriminants . First is "local 
additivity" , i.e. 

M(A 1 +B,A 2 ,...,A N ) = M{A 1 ,A 2 ,...,A N ) + M{B,A 2 ,...,A N ). 

Second is permutation invariance , i.e . 

M(A 1 ,A 2 ,...,A N ) = M(A t{1) ,A t{2) ,...,A t{n) ),t g S N . 

And the third one is easy formula for the rank one case : 

M(xiyi, ■■■■> xnVn) = det(xiyf + ... + x N yJf), 

where (xi,yi; 1 < i < N) are iV-dimensional complex column- vectors . 
Recall that blocks of pa are defined as 

Aij = E A k e ie ]Al, \<i,j<N. 

l<i<k 

Let us rewrite formula (15) as follows : 

QP(P)=--^ E (-1)^" {,7) M(A t(1)i(t(1) ,...,A tW)(tW ); (41) 

ct,t£Sn 

Using this formula (41) we get the following expression for quantum permanent of bipartite 
density matrix pa using "local" additivity of mixed dicriminant in each matrix component : 

QP( P A) = j ri E E M (^i e n(l)4 2 (l)4'-'^ e nw4 2W 4 iv )- 

ti,...,t N Ti,T2G5jv 

Using rank one formula above and formula(40) , we get that 

E M(A tl e n( i ) eJ 2(1) ^ t 1 ,...,A tjv e 7l(JV) eJ 2(JV) 4 JV ) = \M(A tl , A tN )\ 2 . 
Ti,T 2 eS N 



The last formula gives the following , intermediate , identity : 

1 

M 



QP( P a) = -^ E \M(A tl ,...,A tN )\ 2 . (42) 
t\,...,t N 



What is left is to "collect" in (42) , using invariance of mixed discriminants respect to permu- 
tations , all occurances of M(B r ) (as defined in (39)) , where r = (r±, ...,rfc) 6 Ik,N ■ 
It is easy to see that this number N(ri, r^) of occurances of M(B r ) is equal to the coefficient 
of monomial x^x^ 2 ...x r k k in the polynomial (x\ + ... + Xk) N ■ 
In other words , JV(ri, r^) = ri \ Nl rk \ , which finally gives that 

r)| 2 
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Using formula(39) for coefficients of determinantal poiynomial det(J2i<i<k x iAi) we S e t that 
\\Pa\\1 = E l«n,..,rj 2 n!r 2 !...r fc ! = QP(p A ) 

(ri,...,r fe )e4,jv 

IPutting Part 1 and Part 2 together we get in the next corollary a formula expressing permanents 
of positive semidefinite matrices as squares of G-norms of multilinear polynomials . A particular 
, rank two case , of this formula was ( implicitely ) discovered in [3] . 

Corollary A.l: Consider complex positive semidefinite N x N matrix Q = DD^ , where a 
"factor" D is N x M complex matrix . Define a complex gaussian vector z = D£ , where £ is 
an M -dimensional complex gaussian vector as in theorem 3.8 . 

The following formula provides unbiased nonnegative valued random estimator for Per(Q) : 

Per(Q) = E^ N (\ Zl \ 2 ...\z N \ 2 ). (43) 

Proof: Consider the following m-tuple of complex N x N matrices : 

Diag= (Diag 1 ,...,Diag m ))Diag j = Diag(D(l,j), ...,D(N,j)),l <j< M. 
Then Pr>i as (xi, ...,x m ) = IIi<i<7v {Dx)i , where (Dx)i is ith component of vector Dx . 
Thus Part 1 of theorem 3.8 gives that H-FfoiagHo = -E^ 1i ... iZjv (|zi| 2 ...|ztv| 2 ). 

It is easy to see that the block representation of bipartite density matrix poiag associated with 
m-tuple Diag is as follows : 





A lt2 ■ 


■ A 1>N \ 




A 2 ,i 


A 2 ,2 ■ 


■ A 2) N 


■>Ai,j = Q(h j) e i e j 


V Mi 


An,2 ■ 


■ An^n J 





Therefore QP{pr>i a g) = Per(Q). 
Now Part 2 of theorem 3.8 gives that 

Per(Q) = QP(pDiag) = H^biagllc = E Zl) ... )ZN (\ Zl \ 2 ...\z N \ 2 ). (44) 

I 

Remark A. 2: Corollary (A.l) together with a remarkable supermultiplicative inequality for 
the G-norm ( jSJ, ^) give a completely new look at many nontrivial permanental inequalities 
, such as famous Leib's inequality etc, and allow new corellational inequalities for analytic 
functions of complex gaussian vectors and new ( "short" ) characterizations of independence of 
analytic functions of complex gaussian vectors . More on this will be described in ^8 . I 
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B Wick formula 



In the next theorem we recall famous Wick formula (see , for instance , |28] ) . 

Theorem B.l: Consider complex 2N x M matrix A and a real M -dimensional gaussian vector 
x with zero mean and covariance matrix E{xx T ) = I . Define (yi, •••,2/2A r ) T = Ax . Then the 
following Wick formula holds 

W(A)=:E( J] Vi ) = Haf(AA T ), (45) 

l<i<27V 

where hafnian Haf(B) of 2N x 2N matrix B is defined as follows : 

Haf(B)= J2 II B(pu qi) (46) 

i<pi<p 2 <...pjv;pi<giv:Piv<gjv<2A'' i<i<N 



Let us show how formula (43) follows from (45) . 



Proposition B.2: Suppose that complex N x M matrix D in Theorem 1.4 can be written as 
D = C + iB . Consider the following complex 2N x 2M matrix A : 

C + iB iC-B \ 
C -iB -B-iC J ' 




Then 



AA T 



DD^ 
DD^ 



and W(A) = E Zlr „ tZN (\z\\ 2 ...\zjy\ ) , where the expression 

Z N ( 



^^...^jvd-Zil 2 ...^^! 2 ) is the same as in Corollary (A.l) 



As it easy to see that Haf(AA T ) = Per(DD^) , thus , using Wick formula (45) we reprove 
formula (43) . 

Summarizing , we can say at this point that formula(43) is essentially a different way to 
write Wick formula . (We thank A.Barvinok for pointing at this observation and reference 28 ). 
From the other hand formula(43) is a direct corrolary of formula(24) for the case of tuples of 
diagonal matrices . It is easy to see that one can also consider upper triangular matrices . 
More generally consider the following group action on tuples of square complex matrices A = 
(A u ...,A k ): 

A x ,y = {XAiY, ...,XA k Y), det(Xr) = det(X) det(Y) = 1 (47) 

As 

Pa x<y {xx, ...,x k ) = det( XiXAiY) = det(A) det(Y)P A (x 1 , ...,x k ) 

l<i<k 
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this group action does not change corresponding determinantal polynomial . 

Finally , it follows that Wick formula is a particular case of formula(24) when there exist two 

matrices X and Y such that det(XY) = det(X)det(Y) = 1 and matrices (XAiY,...,XA k Y) 

are all upper triangular , or , in Lie-algebraic language , there exists two nonsingular matrices 

X and Y such that the Lie algebra generated by (XAiY, ...,XAf~Y) is solvable . 

It seems reasonable to predict that formula(24) might be of good use in combinatorics desribed 

in HH1 ■ 

C Short probabilistic proof of formula (24) 

C.l Hilbert space of analytical functions 

Consider a Hilbert space L^^g of analytic functions 

f(x\,X2,---,Xk)= J]] a rx,--i r k X l X 2 ■■■ x k ' 
(rx,—,rk) 

where the G- inner product is defined as 

<f,9>G= ari,...,r k b ri ,...,r k riL.r k \ (48) 

It is easy to see that L^g is a closed proper subspace of L2(C k ,/j,) , where \x is a gaussian 
measure on C k , i.e. its density function 

p(z) = i-e-N 2 

Proposition C.l: Suppose that f,g £ Li2(C and the matrix U : C k — > C k is unitary , i.e. 
UU* = I . Then 

< f(Ux),g >£ 2 (cfc, M )=< f,g(U*x) >l 2 (c"s^) 

Proof: This is just a reformulation of a well known obvious fact that p(z) = piUz) (e— \z\ 2 = 
e—\Uz\ 2 ) for unitary U . I 

Lemma C.2: Let P(xi, X2, x k ) be a homogeneous polynomial of total degree N and g £ 
L 2 (C k ,n) . Then for any matrix A : C k — > C k the following identity holds : 

< P{Ax),g > L2(CM = < /, g(A*x) >l 2 (C>°,») (49) 

Proof: First , there is an unique decomposition g = Q + 5 , where Q(x±, X2, x k ) is a 
homogeneous polynomial of total degree N and < R,5 >L 2 (C k ,^) = for any homogeneous 
polynomial R of total degree N . 

As P(Ax) is a homogeneous polynomial of total degree iV for all A thus < P(Ax), S > L 2 (C k ,n)= 
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It is left to prove (49) only when g is a homogeneous polynomial of total degree N . We already 
know that (49) holds for unitary A . Also , because of formula(23) , in this homogeneous 
case (49) holds for diagonal A . To finish the proof , we use the singular value decomposition 
A = VDiagU , where U, V are unitary and Diag is a diagonal matrix with a nonnegative entries 
I 

Remark C.3: A homogeneous part of Lemma has been proved in [1] using the fact that the 
linear space of homogeneous polynomials of total degree N is spanned by N powers of linear 
forms . | 

C.2 Unbiased estimators for Quantum permanents 

Remark C.4: Consider a four-dimensional tensor p(i\, 12, 23, 14), 1 < ii,*2>*3>*4 — N ■ One 
can view it as a block matrix as in (3) , where the blocks are defined by 

A h,h ='■ {p(h,i2,ji,32) ■ 1 < h,h < N},1 < h,ji < N 

We also can permute indices : p(i n m , \(2) > V(3) 1 V(4) ) j an d get another block matrix . The 
main point is that it follows from formula(16) that a permutation of indices does not change 
the quantum permanent QP(p) . In what follows below we will use the following simple and 
natural trick : permute indices and use mixed discriminants based equivalent formula(15) for 
QP(p) based on the corresponding block structure . 
I 

The next proposition follows directly from the definition . 
Proposition C.5: 

1. Consider a block matrix p as in (3) and associate with it the following operator T : 
M(N) -> M(N) , T(X) = Ei<ij<NX(i,j)Ai d . 

Let X be a random complex zero mean matrix such that E(\X(i, j)\ 2 = 1 and for any two 
permutations T\,T2 the set of entries 

{Xij : j = Ti(i) or j = T2(i)} consists of independent random variables . 
Then 

QP(p) = £(det(T(X))det(X)). (50) 

2. Consider a N x N matrix A and a random zero mean vector z £ C N such that E(z{Zj) = 
for all i ^ j ■ Then 

Per(A) = E( \{ (Tz\ ]J z~) (51) 

l<i<N l<i<N 

Let us present now a promised short probabilistic proof of (24 ) : 

Consider without loss of generality a iV 2 -tuple of N x N complex matrices A = (Ann, Aj-jyjv))- 
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Recall that the (BUDM ) associated with the /c-tuple 
A = (Ai, ...,Ak) is defined as : 

PA(n,i2,ji, h) =■ M i i,i2)Ai(j 1 ,j 2 ) 

l<l<k 

Define the following operator T(X) = J2i<ij<N j)Bij , where N x N matrix Bij = 
{pA(i,j,m,l) : 1 < m,l < N} . It is easy to" see t hat T(X ) = C*C(X) , where C(X) = 
J2i<i,j<N j)A(ij) . Thus QP(pa) = Ex(det(T(X)det(X)) , where random gaussian matrix 
X has the density 

p(i)4^ (xr) ' 

I.e. the entries X(i,j) are IID canonical complex gaussian random variables . 
Finally , we get that (using for the first identity (50) and for the second (49) ) 

QP(pa) = E x (det(T(Xjtet(Xj) = E *(\ det(C(X)| 2 ) 
But , E x (\det(C(X)\ 2 ) = \\P A \\ 2 G from (23). 

Similarly , the permanental formula (43) can be proved using (51) . 
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